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NOTATION 


C,.  Frictional  drag  coefficient 

Cp  Pressure  coefficient,  Ap/(  1 /2  p U^) 

Cp  Residual  drag  ccx-fficient 

Total  drag  coefficient 
Prismatic  coefficient  of  forebody 
O Maximum  diameter  of  model 

D|  Flat-faced  diameter  or  friction  drag 

Djj  Residual  drag 

Total  drag 

k Turbulence-stimulator  height 

kjj  Nondimensional  rate  of  change  of  curvature  on  forebody  at  juncture  of  flat  face 

with  transition  curve 

k,  Nondimensional  rate  of  change  of  curvature  on  forebody  at  juncture  with 

parallel  middle  body 

L Length  of  model  or  prototype 

L,  U'ngth  of  fore  body 

L|^  Length  of  afterbody 

R Reynolds  number* 

Ry  Reynolds  number,  R^  = MOIv 

r Nondimensional  radius  of  curvature  at  forward  point  of  body 

S Surface  area  of  model 

Sp  Surface  area  of  forebody 

S,  Surface  area  between  x = x„  and  x = L 

Sjj  Surface  area  of  afterbody 

S^j  Surface  area  between  x = 0 and  x = Xg 

*Uniets  otherwise  defined  all  Reynolds  numbers  are  based  on  the  freestream  velocity  Uq  and  a length 
denoted  by  the  subscript. 
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1 

Sxg-Xo 

Surface  area  between  x = and  x = Xj 

4 

s 

Arc  length  along  a meridian 

Arc  length  of  model  along  a meridian 

* 

T 

Forcbody  slenderness  parameter  T = (Lp/D)/Cp^ 

U 

Potential  flow  velocity  on  the  model 

Uo 

Speed  of  model  or  prototype 

»k 

Laminar  boundary  layer  velocity  at  trip  height  k 

X 

Axial  distance  from  model  forepoint 

Xj 

Length  of  laminar  flow  on  model  forebody 

y 

Forebody  radius 

a 

Tangent  angle  of  model  profile 

• 

Stimulator  drag  coefficient,  AD,/(l/2  p Uj^  2 7ry|^  k) 

ACt 

Stimulator  drag  coefficient,  AD,/(  1 /2  p S)  ^ 

AD, 

Stimulator  drag 

Ap 

Pressure  change  on  model  from  potential  flow 

b* 

Boundary  layer  displacement  thickness 

1 

ri 

Nondimensional  forebody  radius 

• 

0 

Boundary  layer  momentum  thickness 

^xg 

Boundary  layer  momentum  thickness  at  x = Xj 

1 

^XJ-Xo 

d at  X = xg  from  turbulent  flow  originating  at  x = Xj, 

i X 

dU 

Pressure  gradient  parameter,  — — 

1 

1 

Kinematic  viscosity  of  fluid 

1 

1 

1 i 

Nondimensional  axial  distance  from. forward  point 

4 

1 

i 

1 

' p 

Fluid  density 

« 

I 

3 

I 

i 

Subscripts 

j 

• * 

k 

Location  of  stimulator 

i 

, « 

Laminar  flow 

1 o 

Location  of  virtual  origin 

* 
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Turbulent  flow 
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ABSTRACT 


Resistance  has  been  measured  of  nine  bodies  of  revolution,  having  equal 
volume  but  varying  forebody  shapes.  Forebody  shapes  ranged  from  extremely 
blunt  to  extremely  fine  and  included  two  that  were  flat  faced.  The  forebodies 
were  altered  by  changing  their  length-to-diameter  ratios  (L/D’s)  and  prismatic 
cixMficients.  Drag  results  indicate  that  when  the  prismatic  coefficient  is  fixed 
and  the  L/D  is  decreased,  the  residual  resistance  will  increase  modestly. 
Increasing  the  prismatic  coefficient  at  small  L/D’s  increases  residual  resistance; 
however,  at  moderate  L/D’s  it  does  not.  The  results  suggest  that  a fiat-faced 
shape  in  itself  does  not  increase  resistance. 

In  addition  to  resistance  experiments,  transition  regions  on  the  models 
were  located,  using  hot  film  probes.  Calculations  predicted  laminar  separation 
on  five  of  the  model  forebodies.  The  hot  film  measurements  confirmed  that 
separation  did  occur  at  the  locations  predicted;  downstream  of  the  separation 
locations,  turbulent  fiow  occurred  immediately.  The  remaining  four  fore- 
bodies exhibited  well-defined  natural  transition  regions.  Flow  properties  in 
the  transition  regions,  measured  by  the  hot  film  gages  have  been  compared 
with  predicted  spatial  amplification  ratios  of  disturbances,  calculated  by 
linear  stability  theory.  Results  have  failed  to  show  a single  relationship 
between  measured  fiow  properties  and  computed  spatial  amplification  ratios. 
Correlation  of  amplification  factors  with  fiow  regimes  varied,  both  with 
forebody  shape  and  Reynolds  number. 


ADMINISTRATIVE  INFORMATION 

This  work  was  funded  and  authorized  by  the  Naval  .Sea  Systems  Command  under  the 
hydrodynamics  part  of  the  very  high  speed  submarines  program.  Task  Area  SF  434  215  IZ, 
Work  Unit  1-1  520-004.  and  the  improvements  in  submarine  hydrodynamics  program.  Task 
Area  SSL66007,  Work  Unit  1-1552-135. 


INTRODUCTION 

Experiments  have  been  conducted  on  nine  bodies  of  revolution  having  different  forebody 
shapes.  Ihe  primary  purpose  of  the  experiments  was  to  determine  what  effect  changing 
' forebody  shape  would  have  on  resistance.  As  a means  to  this  end,  the  roles  of  transition, 
laminar  separation,  and  turbulence  stimulation  in  the  analysis  of  model  data  have  also  been 
investigated.  These  phenomena  have  been  examined  using  hot  film  probes,  located  in  the 
model  surfaces  - resulting  in  a'  new  method  for  analyzing  model  drag  data  as  reported  by 
McCarthy,  Power,  and  Huang.*  Results  of  both  the  resistance  and  hot  film  measurements  are 
presented  in  this  report. 
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The  experiments  have  been  made  in  two  series.  In  the  first,  the  effect  of  forebody 
length-tcvdiameter  ratio  (L/D)  was  investigated  by  measuring  the  resistance  and  transition 
hK'ations  of  four  models  with  forebodies  having  different  bow  entrance  L/D’s  and  a fixed 
bow  prismatic  coefficient  of  0.667.  The  data  obtained  in  this  series  were  also  used  in 
developing  the  new  method  of  analysis  given  in  Reference  1.  For  completeness,  most  of  the 
results  from  the  Series  1 experiments,  reported  in  Reference  1,  are  repeated. 

In  Series  2.  the  effect  of  very  large  values  of  bow  prismatic  coefficient  on  resistance  was 
investigated.  Properly  shaped  blunt  forebodies  could  create  strong  favorable  pressure  gradients 
in  the  forward  regions  of  the  forebodies  where  transition  could  be  expected  to  take  place  at 
full-scale  Reynolds  numbers.  These  pressure  gradients  would  extend  laminar  How  further  aft 
when  compared  with  flow  on  more  conventional  forebodies,  possibly  producing  a more 
favorable  operating  environment  for  sensors,  in  the  bow  region.  In  Series  2,  five  blunt  fore- 
bodies were  investigated,  including  two  that  were  flat  faced. 

MODELS 

lo  conduct  the  resistance  experiments,  nine  models  were  constructed,  all  bodies  of 
revolution.  Fach  model  consisted  of  a forebody,  parallel  middle  body,  and  fixed  afterbody; 
see  Figure  1.  Seven  of  the  models  had  rounded  forebodies,  and  two  had  flat  faces.  The 
volume  and  the  maximum  diameter  of  ail  models  were  kept  constant.  In  this  section,  the 
method  used  to  generate  forebody  shapes,  the  selection  of  the  forebodies,  and  the  constniction 
of  the  models  will  be  discussed. 

GENERATION  OF  FOREBODIES 
RourKled  Forebodies 

Six  of  the  rounded  forebody  shapes  were  developed  using  the  Granville  family  of 
quadratic  polynomials^;  the  seventh  "rounded”  forebody  was  hemispherical  in  shape.  The 
nondimensional  offsets  of  the  polynomials  and  the  hemispherical  shape  are  given  by  the 
following  two  equations 


’ McCarthy,  J.H.  el  al.,  "The  Roles  of  Transition,  l.aminar  Se/kirarion,  and  Turbulence  Stimulation  in  the 
Analysis  of  A xisvmmetric  Body  Drag,"  Eleventh  Symposium  on  Naval  Hydrodynamics,  London 

^Granville,  P.S.,  “Geometrical  Characteristics  of  Noses  and  Tails  for  Parallel  Middle  Bodies, " NSRDC 
Report  3763  (Dec  1972). 
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Polynomials 


T7  = rK(^)  + k,  K,({)  + 0(£) 


(1) 


where  R(J)  = 2 - I )' 

= 1)^ 

Q({)  = 1 -(^-  !)■*  (4{  + 1) 


Hemisphere 


r?  = 


1/2 


(2) 


where 


wliere 


^ = x/L|  and  tj  = 2y/IJ 

D is  forebody  maximum  diameter 

L|  is  forebody  length 

X is  axial  distance  from  forward  point 

y is  fore  body  radius 

r is  radius  of  curvature  at  ^ = 0 

k|  is  rate  of  change  of  curvature  at  ^ = 1 


By  selecting  proper  values  of  r and  k, , various  parent  shapes  can  lie  derived.  The  values 
of  Lj^  and  D then  determine  the  shape  of  the  particular  forebody.  Once  values  of  r and  k, 
have  been  selected,  the  prismatic  coefficient  ( is  uniquely  determined  and  is  given  by 


r 

15 


(3) 


At  X = U . the  slope  and  the  curvature  ( are  set  at  zero,  matching 

VdxV^^L,, 

their  corresponding  values  on  the  parallel  middle  body.  Ihe  hemispherical  shape  is  fixed 
with  Lj.  = D/2.  At  x = ^ the  slope  matches  that  of  the  parallel  middle  body;  however,  the 
curvature  does  not. 


Flat-Faced  Forebodies 


I he  curves  connecting  the  Hat  lace  to  the  parallel  middle  body  of  the  two  flat-faced 
shapes  were  developed,  using  both  an  elliptical  contour  and  a contour  derived  from  the 
Granville  two-parameter  cubic  polynomial.^  I he  nondimensional  elliptical  contour  offsets  are 
given  by  Fquation  (2),  and  the  Granville  cubic  polynomial  offsets  are  given  by 


K„({)  = 1)^ 

0(£)  = 1 (4f  + 1 ) 


where  now 


{ = X L|  and  r? 


2y  -I), 

D - D, 


/d^x\ 

\dyV,= 


(4) 


2L, 


D-D 


X = L. 


where  Uj  is  the  diameter  of  the  Hat  face. 

l or  the  CJranville  shape,  the  slope  and  curvature  of  the  connecting  curve  match  the  slope  and 
curvature  of  the  Hat  face  and  parallel  middle  body  where  the  curves  join,  and  the  shape  of  a 
parent  body  is  determined  by  values  assigned  to  and  k, . Ihe  values  of  1.,  and  (H-  Dj 4/2 

then  determine  the  contour  of  the  particular  forebody.  The  elliptical  transition  curve  (non- 
dimensional)  is  fixed,  and  variation  in  model  forebody  shape  is  a function  of  only  L|  and 
(D-  D|  )/2.  Ihe  slopes  match  at  the  flat  face  and  parallel  middle  body;  however,  the 
curvatures  do  not. 


/ 

^Granville,  P.S.,  "Geometrical  Characteristics  of  Flat-FaccJ  Bodies  of  Revolution.  " NSRDC  Report  37 1 0 
(Nov  1971). 
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SELECTION  OF  SHAPES 

HDr  the  first  scries  ol  experiments,  lour  models  were  selected.  Since  the  purpose  of  this 
series  was  to  investigate  the  effects  of  I.,  /I),  only  this  parameter  was  varied.  All  the  models 
had  rounded  forehodies  with  prismatic  c(H;fficients  et|ual  to  0.667.  Ihe  first  forebody  model 
selected  (4620-.^)  is  designated  as  the  parent  model,  Lj  /I)  = I.S2;  the  resistance  of  this  model 
is  then  used  as  a i)aseline  to  assess  the  merits  of  the  other  forehodies.  Ihe  remaining  three 
forehodies  have  l.j..  /D  ratios  of  0.5,  1.0,  and  3.0.  The  I.j.  /D  = 0. 5 forebody  has  a hemispher- 
ical shajic'.  Ihe  other  three  have  Ciranville  shapes  and  were  derived  from  a common  parent, 
i.e.,  r = O.X333  and  kj  = 10.0.  Figure  2 shows  the  contours  of  the  four  forebody  shapes. 

Of  the  live  forehodies  selected  for  experiment  in  Series  2,  three  had  rounded  shapes. 
They  were  selected  during  a parametric  investigation  of  forebody  shape,  using  calculated 
pressure  distributions  and  transition  locations  for  guidance.  Pressures  were  computed  for 
several  model  shapes,  using  the  Hess-Smith,  potential-llow-computer  program.'*  Boundary- 
layer  stability  calculations  were  made  for  a few  selected  shapes  at  typical  full-scale  Reynolds 
numlxTS,  using  the  Smith-Gamberoni  method.*-^  Results  were  then  compared  with  similar 
calculations  for  the  parent  forebody  shape.  Model  46  20-.T  .Shapes  selected  were  chosen  for 
their  persistence  of  laminar  How  in  the  full-scale  Reynolds  number  range  I hey  were  derived 
from  a common  parent  - r = 3. 167.  Cp|.  = 0.850.  kj  = 5.0  - having  l.j.  /I)  ratios  of  0.50,  1.00, 
and  1.82. 

The  lirst  llal-faced  forebody  selected  has  a relatively  small  prismatic  coefficient,  0.823. 
and  a forebody  Lj  /I)  of  1.216.  The  connecting  curve  between  the  Hat  face  and  the  parallel 
middle  body  on  this  model  has  an  elliptical  contour.  In  contrast,  the  second  llat-faced  fore- 
body is  extremely  blunt  with  a prismatic  cix'fficient  eipial  to  0.633.  \ (iranville  twev 
parameter  cubic  polynomial  ( k^,  = I.O,  k,  = 3.0)  was  used  for  the  transition  curve  on  this 
model,  fhe  contours  of  the  five  forehodies  of  the  second  series  are  shown  in  Figure  3.  and 
tile  particulars  of  the  nine  mtidels  selected  are  given  in  Fable  I. 

Incipient  cavitation  speeds  have  been  calculated  for  the  nine  forebody  model  shapes;  all 
are  listed  in  Fable  2.  Ihe  model  with  the  earliest  predicted  cavitation  inception  is  4620-7. 


^Iless,  J.I..  and  A.M.O.  Smith,  '‘Calailaiion  of  Potential  Flow  about  Arbitrary  Bodies,"  Progress  in 
Aeronautical  Sciences,  Vol.  8,  Pergairion  Press,  Inc.,  New  York  (1666). 

*Sinith,  A.M.O.  and  N.  (lamberoni,  "Transition  Pressure  Gradient  and  Stability  Theory,"  Douglas  Aircraft 
Company  Report  r.S-26.)88  (Aug  1656). 

^Wa/./jtn,  A.R.  et  al.,  "Spetml  and  Temporal  Stability  Charts  for  the  Flakner-Skan  Bourulary-I.ayer  Profiles, 
Douglas  Aircraft  romjrany.  Report  DAC-67086  (55ep  1668). 


5 


TABLE  1 - MODEL  PARTICULARS 


Model 

B 

D 

D 

L 

D 

_S_ 

£e 

SERIES  1 

4620  1 

0.667 

0.500 

0 

mgm 

29.54 

1.57 

4620-2 

0667 

1.000 

0 

■B 

29.60 

2.67 

4620  3 

0.667 

1.820 

0 

11.17 

29.81 

4.61 

46204 

0.667 

3.000 

0 

11.57 

30.19 

7.46 

SERIES  2 

4620  5 

0.850 

1 000 

0 

10.72 

29.47 

3.21 

4620-6 

0 

10.84 

29.53 

5.46 

4620  7 

■■■ 

0 

10.64 

29.48 

1.90 

4620  8 

■■■ 

0.540 

10.78 

29.63 

3.81 

4620  9 

1.000 

0.500 

10.63 

29.43 

3.35 

D“  0.624  m.  Lp/D«4.400,  Sp/D^  = 40.21;  Volume/D^  » 6.786. 


TABLE  2 - COMPUTED  INCIPIENT  CAVITATION 
SPEEDS  OF  NINE  FOREBODIES 


Incipient  Cavitation  Speeds  in  Meters  per  Second 
for  Water  Depths  of  - 


"mm 

15m 

30  m 

60  m 

90  m 

150  m 

1 

-0  771 

25 

32 

50 

50 

64 

2 

-0  402 

35 

44 

58 

70 

89 

3 

-0  185 

51 

65 

86 

103 

130 

4 

-0  094 

72 

91 

121 

144 

183 

5 

-0  778 

25 

32 

42 

50 

63 

6 

-0428 

34 

43 

57 

67 

85 

7 

- 1 402 

19 

24 

31 

38 

47 

8 

-0.610 

28 

36 

47 

57 

72 

9 1 

-0988 

22 

28 

37 

44 

56 

Note:  It  .'ssumed  that  forebodies  have  negligible  roughness. 


At  a depth  of  15  m this  shape  had  an  inception  speed  of  18.5  m/s.  The  cavitation  calculations 
were  based  on  the  computed,  potential  How,  pressure  distributions  for  smooth  forebodies. 
Rough  forebodies  would  have  lower  inception  speeds.  The  pressure  distributions  on  the 
models  of  the  first  series  arc  given  in  Figure  4 and  of  the  second  series  in  Figure  5.  Minimum 
values  of  the  pressure  distributions  are  listed  in  Table  2. 

CONSTRUCTION 

The  models  were  constructed  of  molded  Fiberglas  and  consisted  of  three  sections.  The 
second  and  third  sections  were  common  to  all  models  and  consisted  of  part  of  the  parallel 
middle  body  and  the  afterbody,  respectively.  The  first  section  included  the  forcbody  and  the 
remainder  of  the  parallel  middle  body.  For  each  of  the  first  sections,  the  length  of  the 
included  parallel  middle  body  was  adjusted  to  insure  constant  volume  for  all  models.  The 
third  section  was  obtained  from  an  existing  model  of  the  research  submarine  ALBAfORE 
(AGSS-569)  by  cutting  the  model  in  two  parts  at  the  position  of  maximum  diameter. 

Special  care  was  taken  in  the  construction  of  the  forebodies  of  the  models.  Profile 
tolerances  were  held  to  ± 0.4  mm;  all  imperfections  were  removed  and  meridians  were  fared. 
The  Fiberglas  was  polished  to  a 0.64-micron,  root-mean-square,  surface  Finish.  This  insured 
roughness  Reynolds  numbers  two  orders  of  magnitude  lower  than  the  number  considered 
necessary  to  trip  turbulence  prematurely. 

EXPERIMENTS 


GENERAL 

I he  model  drag  experiments  were  conducted  in  the  David  W.  Taylor  Model  Basin,  using 
Carnage  2.  Ibis  towing  basin  is  846  m long.  15.5  m wide,  and  6.7  m deep.  Drag  and  speed 
data  were  collected  for  each  model  in  the  speed  range  between  0.51  and  515  m/s  at  intervals 
of  approximately  14  m/s.  A waiting  time  of  5 min  between  runs  was  selected  because  no 
effects  of  residual  turbulence  from  preceding  runs  could  be  detected  after  this  time.  Flush- 
mounted  hot  films  were  inserted  in  the  model  surface  at  selected  locations  to  determine  the 
extent  of  laminar  and  transitional  flow  and  the  location  of  laminar  separation  when  it 
occurred.  Water  temperatures  were  monitored  throughout  the  test  period  to  obtain  water 
densities  and  kinematic  viscosities. 
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RESISTANCE  MEASUREMENTS 

Total  drag  has  been  measured  using  the  standard  DTNSRDC  block-gage-type  dynamo- 
meter. mounted  internally,  on  each  model  and  connected  to  the  carriage  with  a single  stream- 
lined strut.  This  dynamometer  has  an  accuracy  of  approximately  1 percent  at  its  maximum 
load  capability  of  .‘'.‘i  k.  The  model  centerline-submergence  depth  was  fixed  at  2.74  m below 
the  water  surface  for  all  runs.  The  models  were  ballasted  for  a slight  amount  of  negative 
buoyancy  with  zero  moment  at  the  towing  strut.  Previous  work  has  shown  that  the  towing- 
strut  model,  interference-drag  coefficient  is  less  than  0.01  • 10"^  for  the  size  of  models 
investigated.  Since  this  interference  drag  is  within  the  experimental  accuracy  of  model  drag 
data  reported  here,  no  corrections  have  been  made  for  it. 


TRANSITION  MEASUREMENTS 

Table  3 gives  the  locations  of  the  hot  films  in  the  various  models.  Selection  of  the  hot 
film  locations  was  assisted  by  inspection  of  pressure  distributions,  experience  from  previous 
tests,  and  some  trial  and  error.  When  laminar  separation  was  predicted,  a hot  film  was 
inserted  immediately  fore  and  aft  of  the  predicted  location.  The  Curie  and  Skan  modified 
Thwaites  criterion^  of  X = -0.09  was  used  to  predict  the  location  of  laminar  separation  where 
X is  a pressure  gradient  parameter  given  by 


dU 
V ds 


where  d is  boundary  layer  momentum  thickness 
i>  is  kinematic  viscosity 
U is  potential  flow  velocity  on  the  body 
s is  arc  length  along  a meridian. 

The  values  of  0 were  computed  using  the  Granville  integral  method  for  bodies  of  revolution.* 
I he  computed  positions  of  laminar  separation  are  list;d  in  Table  4. 


^Curle,  N.  and  S.W.  Skan.  “Approximate  Methods  for  Predicting  Separatkm  Properties  of  Laminar  Boundary 
l.ayers,“  Aeronautical  Quarterly,  Vol.  8 (1957). 

*Granville,  P.S.,  "The  Calculation  of  the  Viscous  Drag  of  Bodies  of  Revolution,"  David  Taylor  Model  Basin 
Report  M9  (1953). 
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TABLE  3 - LOCATIONS  OF  HOT  FILM  GAGES 


Location 

Nondimensional  Axial  Distance  from  Model  Forepoint  to  Hot  Film  Gages  - 

x/D 

4620-1 

4620  2 

46203 

46204 

4620-5 

46206 

4620-7 

46208 

4620  9 

1 

0.16 

0.29 

0.29 

0.29 

0.30 

0.41 

0.13 

0.05 

0.06 

2 

0.37 

0.45 

0.45 

0.47 

0.42 

0.52 

0.21 

0.17 

0.14 

3 

0.49 

0.61 

0.61 

0.65 

0.54 

0.62 

0.34 

0.29 

0.27 

4 

0.57 

0.86 

0.81 

1.26 

0.75 

0.72 

0.54 

0.46 

0.39 

5 

0.81 

0.98 

1.02 

1.59 

0.82 

0.66 

0.59 

6 

1.22 

1.42 

1.59 

1.91 

1.02 

7 

1.63 

2.04 

1,71 

2.24 

1.23 

• 

8 

2.78 

2.94 

1.87 

3.64 

9 

2.04 

The*  hot  film  probes  were  made  by  Linptronics  Laboratory.  The  active  elements 
(platinum)  were  approximately  1.6  mm  long,  fixed  to  one  end  of  a Pyrex  cylinder  12.7  mm  in 
length  and  2.4  mm  in  diameter.  The  cylinders  were  carefully  mounted  in  predrilled  holes  in 
the  model  surfaces  and  were  held  in  place  with  rubber  cement.  Any  ridge  between  the  model 
surface  and  the  probes  was  held  to  less  than  0.025  mm,  and  the  holes  in  the  model  were 
drilled  perpendicular  to  the  surface  ±0.50  deg.  As  a further  precaution,  each  probe  was 
Iwated  along  a dillercnt  meridian  streamline,  to  avoid  possible  spurious  signals  which  might 
result  from  tripping  of  turbulence  by  probes,  located  at  successively  further  forward  positions. 

Five  channels  of  instrumentation  were  available  and  interchangeable  between  hot  films, 
allowing  the  five  most  significant  probes  to  be  activated  during  a given  run.  Each  channel 
consisted  of  an  anemometer,  linearizer,  oscilloscope,  and  attentuator.  For  each  speed.  20-s 
samples  of  hot  film  output  were  recorded  on  tape.  The  output  of  each  channel  was  kept 
close  to  but  less  than  the  1.4-V  overload  capacity  of  the  recorder  by  adjusting  the  overheat 
ratios  of  the  films,  the  linearizer  amplifiers,  and  the  attenuators.  Overheat  ratios  were  kept 
between  1.035  and  1.050  for  most  of  the  tests;  however,  for  some  of  the  low-speed  runs, 
slightly  higher  ratios  were  needed.  Since  these  adjustments  had  to  be  made  quickly  as  the 
carriage  moved  down  the  basin,  no  record  of  amplifications  or  overheat  ratios  could  be  kept. 
Thus,  the  information  obtained  from  the  hot  films  was  primarily  qualitative  and  was  used  to 
determine  the  nature  of  the  flow  at  the  hot  film  locations.  Figure  6 shows  one  channel  of 
the  instrumentation. 
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TABLE  4 - LOCATIONS  OF  WIRE  STIMULATORS  AND  TRANSITION 
OR  LAMINAR  SEPARATION  WITHOUT  STIMULATORS 


Model 

Wire 

Diameter 

mm 

Wire 

Locations 

x,/D 

Locations  of 

Laminar  Separation  (LS) 

or 

Transition  (T) 

4620  1 

0.61 

0.107 

LS:  x/D- 0.472 

0.322 

0.61 

0.218 

LS:  x/D  =0.894 

0.382 

0.545 

4620-3 

0.25 

0.223 

T: 

0.559 

0.61 

0.223 

1.56  < x/D  < 2.13, 

0.390 

0.559 

33.6  > Rp  • 10‘®>  3.1. 

1.02 

0.223 

0.559 

46204 

0.61 

0.231,0.405.  0.579 

T: 

0.405.  0.579 

1 .26  < x/D  < 3.39, 

0.579 

33.6  > Rp  x/D  10-®  > 3.1. 

4620  5 

0.61 

0.311 

LS:  x/D  =0.371 

46206 

0.61 

0.542 

T:  0.57  < x/D  < 1.04, 

30.6  > Rp  • 10*®  > 3.1. 

4620  7 

0.61 

0.160 

LS:  x/D -0.184 

4620^8 

0.61 

0.078 

T:  0.07  < x/D  < 0.70, 

1 

30.6  > Rq  • 10*®  > 3.1 

LS:  x/D  - 1.25 

46209 

0.61 

0.181 

LS:  x/D -0.202 
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As  shown  in  Figure  7,  analysis  of  the  hot  film,  output  signals  revealed  that  four  different 
types  of  boundary-layer  flow  could  be  distinguished: 

1.  Laminar  flow,  when  hot  film  signals  were  of  steady  amplitude. 

2.  A smooth  wavelike  disturbance  that  had  identifiable  frequency  content,  superimposed 
on  a laminar  flow,  and  no  doubt  associated  with  Tollmien-Schlichting  waves.’  (An  early 
(Type  2a)  and  advanced  stage  (Type  2b)  of  wavelike  How  is  illustrated  in  Figure  7.) 

3.  An  intermittent  turbulent  bursting  flow,  reported  by  Schubauer  and  Klebanoff,'® 
preceded  and  followed  in  time  by  periods  of  Type  (2)  flow. 

4.  A fully  turbulent  flow  of  random  nature. 

Under  fixed  test  conditions,  probe  outputs  indicated  that  Types  (2)  or  (3)  and  Types  (3)  or 
(4)  could  alternately  occur  on  repeat-run,  tape-recorded,  time  samples  of  the  signal  at  a given 
probe  location.  For  a given  location  of  a hot  film  probe  on  a forebody,  the  speed  at  which 
a particular  type  of  flow  occurred  could  be  determined  to  within  1/4  m/s. 

TURBULENCE  STIMULATORS 

Resistance  of  each  model  was  measured  both  with  and  without  turbulence  stimulators. 

In  the  Series  1 experiments.  Model  4620-3  was  used  to  evaluate  stimulators  of  the  following 
types,  sand  strips,  wires,  and  studs.  The  results  of  these  experiments  have  been  reported  in 
Reference  1 and  will  not  be  discussed  here.  The  0.61-mm  wire  was  chosen  for  all  of  the 
remaining  resistance  tests  because  it  proved  to  be  an  effective  trip  in  the  present  experimental 
arrangement,  having  a precisely  defined  geometry  consistent  from  test  to  test  and  being  easy 
to  install  and  remove.  More  parasitic  drag  data  were  also  available  for  the  wire  than  for  the 
other  types  of  trips. 

Table  4 gives  the  wire  sizes  used  and  the  locations  at  which  the  wires  were  installed  on 
the  nine  models.  With  the  exception  of  Model  4620-3,  only  0.61-mm  diameter  wires  were 
used.  Whenever  possible,  experiments  were  conducted  with  wires  installed  at  x/L  = 0 05,  the 
traditional  location  used  at  the  Center  for  trip  wires,  where  L is  the  model  length.  Exceptions 
(Kcurred  when  laminar  separation,  which  acts  as  a natural  trip,  or  transition  took  place  forward 
of  this  location.  When  laminar  separation  occurred  forward  of  x/L  = 0.05,  a trip  was  placed  a 
short  distance  ahead  of  the  separation  location. 

’Schlicting,  H.,  "Boundary  Layer  Theory,"  McGraw-Hill,  Inc.,  New  York,  Chapters  12,  16,  and  17  (1955). 

'“Schubauer,  G.B.  and  P.S.  Klebanoff,  "Contributions  on  the  Mechanics  of  Boundary  Layer  Transition." 

National  Advisory  Committee  for  Aeronautics  TN  3849  (1955). 
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REDUCTION  OF  RESISTANCE  DATA 


I he  value  of  the  residual  drag  coefficient  / 1 /2  p S is  used  to  determine  the 

relative  merits  ot  different  models,  where  Dj^  is  the  residual  drag;  p,  the  density  of  the  fluid; 
the  speed  of  the  model;  S,  the  model  surface  area.  In  the  present  experiments,  values  of 
have  been  determined,  using  a method  developed  in  Reference  1.  In  contrast  to  traditional 
methods,  the  new  method  takes  into  account  the  location  of  transition  on  the  model  in  the 
analysis.  This  is  important  when,  as  in  the  present  case,  the  object  is  to  evaluate  competing 
models  with  vastly  different  forebody  shapes  and  values  with  relatively  small  differences. 

I he  extent  of  laminar  flow  can  vary  greatly  on  such  models  as  the  forebodies  vary  from  blunt 
to  tine.  This  inlluences  the  measured  model  resistance  and  if  not  recognized  will  significantly 
.liter  the  relative  body  drags.  An  outline  of  the  method  is  given  as  follows. 

In  the  analysis  it  is  assumed  that  two  types  of  flow  regimes  exist  on  a model  hull 
laminar  ^ and  turbulent ,.  When  a significant  length  of  transitional  flow  occurrs,  the 
longitudinal  extent  xj  of  the  laminar  flow  regime,  having  a wetted  area  will  be  assumed 
to  terminate  where  intermittent  turbulent  bursting  first  occurs;  see  Type  3 of  Figure  7.  When 
the  wire  stimulator  is  installed,  it  is  assumed  that  xg  will  coincide  with  the  wire  location. 

When  laminar  separation  occurs,  it  is  assumed  that  Xg  coincides  with  the  computed  separation 
location.  The  latter  two  assumptions  are  in  accord  with  available  experimental  data.' 

ITe  value  of  C’^  can  be  expressed  as 

Cp  = C.^-Cp- AC.^  (6) 

where  C'.,  is  the  total  drag  coefficient  DjjMlp  Uj  S 

C,  is  the  frictional  drag  coefficient  l>yj\l2  p\J^  S 
('.j.  is  the  wire  drag  coefficient  AD.j,/l/2  p S 

and  D.J.  is  the  total  drag 

D,  is  the  frictional  drag 
AD.p  is  the  stimulator  drag 

S is  the  model  surface  area 

Ilie  value  of  C'.j^  is  obtained  by  measurement;  however,  the  value  of  C,,.  and  AC.j.  must  be 
obtained  from  empirical  relations. 

Following  Reference  1,  the  drag  coefficient  ACp  of  the  wire  has  been  taken  as  0.75  for 
all  experimental  conditions  where  AC^  is  defined  as 
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(7) 


ACy  = AD^/ (ip  11^2  27ry,^k) 

where  k is  the  wire  diameter 

U|^  is  the  velocity  in  the  laminar  boundary  layer  at  height  k 
yi^  is  the  radius  of  the  model  at  the  wire  location 

rtie  relationship  between  AC^j  and  AC.j.  is 

lu^ J — ^coso,^ 


(8) 


where  Q|^  is  the  tangent  angle  of  the  model  profile  at  the  wire  location.  Values  of  U|^  were 
determined  from  the  Pohlhausen  one-parameter  family  of  velocity  profiles,  using  laminar 
boundary-layer  parameters,  calculated  by  the  method  given  in  Reference  8. 

To  calculate  Cj  . a virtual  origin  < xg  is  defined  so  that  the  frictional  drag  of  a 
turbulent  How  over  the  wetted-surface  area  of  the  model  between  x = x^  and  x = L is  equal 
to  the  total  frictional  drag  on  the  model  (laminar  and  turbulent)  plus  the  drag  of  a stimulator, 
if  present.  This  situation  is  illustrated  in  Figure  8 and  can  be  expressed  by 


D,  +ADr  = D,.^(SL.,^) 


(9) 


Kquation  (6)  can  now  Ir,-  written  as 


where 


I (S|  -X  ^ 

I ( L - X„ 


(10) 


By  analogy  with  bquation  (9)  and  by  converting  to  coefficient  form,  we  can  also  write 


S 


^1,  <^xti-x„> 


xg-x„ 


5 + AC..J. 


(11) 


This  expression  provides  a means  to  evaluate  x^y 

Ihe  frictional  drag  arising  from  the  turbulent  flow  region  on  the  model  is  now 

considered  to  be  the  drag  due  to  the  flow  over  an  equivalent  flat  plate  of  constant  width  - 
2jry  “ S|  x„).  With  this  assumption,  frictional  drag  becomes  a function  of  Reynolds 

number  alone,  and  Equations  ( 10)  and  (11)  reduce  to 
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(12) 


.i!id 


’L  - X, 


D 

*^Xg-Xo 

^'l  ,<*^xe-Xo>  "17  = p 


(13) 


\\licre  0 = ( . Xy /K|  )/(S^g  _ yS).  Also,  using  flat-plate  theory,  we  can  write  for  the  flow 
o\er  the  area  „ 


l ,<^xg-x„' 


(14) 


combining  hquutions  ( 13)  and  ( 14)  we  get 


*^xp 

“s~  ^ 


(15) 


I he  empirical  relation  used  to  evaluate  the  turbulent  frictional  drag  is  the  Schoenherr 
Oat-plate  formulation" 

v/c^  log, 0 (R^  C,,)  = 0.242  (16) 


where  R^^  is  a Reynolds  number  based  on  the  length  of  the  plate.  Using  Equations  (16)  and 
(14)  and  changing  to  the  present  notation,  we  can  write 


.Substitution  for  C,  ^ in  Equation  (14)  yields 

K,„  = R,,-.M,,5,  p (XOx.-xj; 


(17) 


(18) 


lo  obtain  x^,.  Equation  (15)  i.s  solved  iteratively,  after  initially  setting  j3=  1.0  to  obtain 
R^  . I'hen  x^  can  be  obtained  directly  from  Equation  (18).  The  value  of  U|,p(Sxp)  in 
Equation  ( 15)  is  obtained  using  laminar  boundary-layer  theory  for  axisymmetric  bodies. 


"5kht)€nherr,  K.E.,  "Reshtance  of  Flat  Plates  Moving  through  a Fluid,"  Tranuctioni  of  Society  of  Naval 
Architects  and  Marine  Engineers,  Vol.  40  (1932). 
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Tlic  expression  is 


Once  the  value  of  Xj,  is  determined,  C'l.^  (S,  ^nd  Cj^  can  be  obtained,  using  Hquations 
( lo)  and  ( 10).  fhc  value  of  obtained  should  be  the  same  for  both  tripped  and  untripped 
( AC^  = 0)  models. 

When  trips  are  located  near  the  leading  edge,  and/or  the  trip  has  a high  parasitic  drag, 
the  computed  value  of  x„  can  become  negative.  When  this  occurs,  the  following  approxima- 
tion is  made  when  computing  C’j^  from  bquation  ( I 2) 

^1.  • X(,  L “ Xj, 

— — ___ 

S ~ L 

Ihe  results  obtained  in  .Series  I experiments  were  computed  using  the  previously 
described  equations.  However,  in  Series  2 experiments,  whenever  x was  directly  involved  in 
the  calculations,  it  was  replaced  by  the  arc  length  along  the  meridian  s.  Thus  the  virtual 
origin  was  located  at  a position  s^,  from  the  leading  edge  of  the  plate  of  width 
2ny  = Sl  “^o*’  operating  at  a Reynolds  number  KsL-s^,  ~ where  s^  is 

the  total  arc  length  of  the  model.  This  change  was  necessary  to  obtain  results  on  the  two 
llat-faced  models.  To  compare  results,  values  were  computed  using  both  s and  x for 
Models  5 through  7.  There  was  a difference  of  approximately  I percent  in  the  values  of  C',^ . 

RESULTS 

TRANSITION  LOCATIONS 

On  Models  I,  2,  5,  7,  and  *),  laminar  separation  was  predicted  to  occur  on  the  forebody 
sections  of  the  models  at  the  locations  given  in  Table  4.  For  the  Reynolds  number  range  used 
in  these  experiments,  hot  films  located  immediately  before  the  predicted  locations  indicated 
laminar  flow.  Hot  film  gages,  located  immediately  aft  of  the  predicted  locations,  indicated  a 
disturbed  or  turbulent  How.  It  was  concluded  from  these  measurements  that  transition  was 
occurring  between  these  two  gages  and  that  separation  was  tripping  turbulent  flow,  thus 
justifying  the  use  of  the  predicted  laminar  separation  location  as  the  transition  location. 


Models  3,  4,  6.  and  8 exhibited  no  laminar  separation.  The  flow  regimes,  measured  by 
the  hot  film  gages,  are  given  in  Figures  9 through  12.  Experimental  data  have  been  compared 
with  predicted  spatial  amplification  ratios  of  disturbance  as  calculated  by  linear  stability 
theory.  These  ratios  were  calculated  using  a computer  program  developed  at  the  Center, 

DABL  (disturbance  amplification  in  boundary  layers).*'^  In  Figures  9 through  1 2.  predictions 
corresponding  to  spatial  amplification  ratios  of  e^,  e*,  c^,  e^,  e” , and  e'^  are  compared  with 
tlow  regimes.  The  results  given  in  these  Figures  fail  to  show  a single  relationship  between  the 
measured  flow  properties  in  the  transition  regions  of  the  four  models  and  the  corresponding 
computed  spatial  amplification  ratios.  Correlations  of  amplification  factors  with  flow  regimes 
vary  both  with  forebody  shape  and  Reynolds  number,  fable  5 gives  the  principal  results 
obtained  for  each  model. 

It  has  been  assumed  in  this  report  that  flow  transition  originates  with  the  start  of 
turbulent  bursting  (Type  3 flow  of  Figure  7).  The  critical  amplification  factor  correlating 
with  this  flow  transition  condition  varies  considerably  from  model  to  model  and  to  a lesser 
extent,  on  individual  models,  with  changes  in  Reynolds  number.  On  Model  3 this  variation  is 
estimated  to  change  from  e'^  at  low  values  of  R^  to  e**  at  the  highest  values  of  R^^.  The 
corresponding  values  on  Model  4 are  from  e**  to  e^;  on  Model  6,  from  e**  to  e^*^.  The  data 
available  for  Model  8 indicate  a nearly  constant  value  of  e’.  However,  due  to  the  failure  of 
a hot  film  probe,  data  at  the  critical  location  at  the  higher  values  of  R,j  are  lacking  for  this 
model. 

RESISTANCE  RESULTS 

The  C'n’s  obtained  for  each  model  are  plotted  in  Figures  13a  through  13i  as  functions  of 
Reynolds  number  Rj . fhe  dotted  line  in  each  figure  represents  the  faired  mean  values  of 
obtained  from  all  of  the  experiments  with  a given  model,  both  with  and  without  turbulence 
stimulators  installed.  The  s<did  lines  roughly  define  /.ones  of  scatter  at  Reynolds  numbers 
greater  than  R^^  = 10^.  The  scatter  ranges  as  much  as  tO.O.SO  • 10*^  about  the  mean  value 
for  some  of  the  models. 

•Stability  results  from  this  program  for  Models  3 and  4 (Figures and  10)  differ  from  stability  results  given 
in  Reference  1 for  the  same  two  models  (Figures  6 and  7 of  that  Reference).  Results  given  in  Reference  1 
were  obtained  by  hand  calculations,  using  charts  from  Reference  6.  The  discrepancy  no  doubt  results  from 
approximations  involved  in  using  interpolations  required  in  determining  the  amplification  factor,  using  hand 
calculations. 

'^von  Kerc/ek,  C.  and  N.C.  Groves,  "Disturhance  AmpUjication  in  Boundary  l.ayen,"  paper  in  preparation. 
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T ABLE  5 - CORRELATION  BETWEEN  AMPLIFICATION 
FACTORS  AND  MEASURED  TRANSITION  DATA 


Amplif  icalton 
Factor 

Model  3 

Model  4 

1 

Model  6 

Model  8 

,'3 

?2  < Ry  • 10  ^ ' 35, 
fully  tu'bulerd 

5<  Rq  • 10  ^ < 22. 
nearly  fully  turbulent 

20<  Rp  • 10  ^ •-  35, 
fully  turbulent 

Probably  fully  turbulent  at 
all  values  of  Rq  tested 

18  < Rp  • 10  ® < 32. 

iTetween  turbulent 

and  intermittent 

5<  Rp  • 10  ® < 18, 
intermiffenr 

5<  Rp  • 10  ® < 32, 
fully  turbulent 

tl 

e 

27<  Rp  • 10  ^ < 35, 
fully  turbulent 

13<  Rq  • 10  ^ <-  27, 
intermtttent 

5<  Rq  • 10'^  < 13, 
wavelike 

20<  Rp  • 10'^  < 35, 
fully  turbulent 

9<  Rp  • 10  ® < 32. 
intermittent 

5<Ro  • 10  ® -'9, 
wavelike 

5 • Rp  • 10  ® < 32, 
fully  turbulent 

20<  Rp  • 10  22, 

intermittent 

5<  Rjj  • 10'^  • 20. 

wavelike 

26<  Rp  ■ 10  ® •'  35, 
fully  turbulent 

18  < Rp  ■ 10  ® < 26, 

jntermiTtenf 

24  < Rp  • 10'®  ^ 32, 

between  inter mittent 

and  wavelike 

5<Ro  • 10  ®<24, 
wavelike 

5<  Rp  • 10  ® < 32, 
fully  turbulent 

Close  to  intermittent 
at  Rq  = 12 

10  < Rjj  • 10  *■  < 20. 
wavelike 

5<  Rq  • 10"^  < 10, 
lamir^ar 

30  < Rp  • 10  ® < 35. 
tnte'^mittent 

15<  Rp  • 10  ® <30, 
wavelike 

5 < Rp  • 10'®  < 15. 
laminar 

5<  Rp  • 10  ® < 32. 
wavelike 

5<  Rp  • 10'®  < 32, 
intermittent 

^0  ’ ^0  ^ ' 10.  wavelike 

Rt,  * 10  10. 

laminar 

Probably  wavelike  in  ranqe 
given  by  10<  Rp  • 10'^  < 35 

10<  Rp  • 10  ® <23, 
laminar 

Probably  Iammar  at  all 
values  of  Rq  tested 

23<  Rp  • 10  ® < 32. 
probably  wavelike 

5<  Rp  • 10  ® <23, 
Iammar 

5<  Rp  • 10  ®<  11, 
wavelike 

e3 

Probably  iammar  at  all 
values  of  Rq  tested 

6<  Rp  • 10  ®<  16, 
Iammar 

Probably  Iammar  at  all 
values  of  Rq  tested 

5 < Rp  • 10  ® < 12, 
Iammar 

Probably  Iammar  at  all 
values  of  Rq  tested 

5<  Rp  • 10  ® < 8. 
Iammar 

Probably  Iammar  at  all 
values  of  Rq  tested 

At  Reynolds  numbers  less  th;in  1.0  • lO’,  the  spread  of  the  data  becomes  much  larger.  This 
results  primarily  from  the  greater  inaccuracy  involved  in  measuring  low  drags  in  the  low-speed 
range.  Figures  14a  through  14j  give  the  values  of  the  measured  both  with  and  without  a 
turbulence  stimulator  installed  for  the  nine  models. 

As  stated  previously,  the  analysis  method  described  in  Reference  1 was  used  to  compute 
the  values  of  C'^ . Ideally.  C'^’s  obtained  from  experiments  with  both  tripped  and  untripped 
• models  should  be  equal  when  this  method  is  used.  However,  an  inspection  of  Figure  13  shows 

that  a small  difference  of  as  much  as  0.03  • 10'^  between  the  tripped  and  untripped  values 
^ . of  C^  can  be  discerned  for  roughly  one-third  to  one-half  of  the  model  data. 

i 
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Wlien  a discrepancy  is  evident,  the  bare-hull  data  give  lower  values  of  than  do  the  tripped 
data.  On  Models  4,  6,  7,  and  ^ for  Rj  > 0.6  • 10^  and  on  Models  .1,  ,S,  and  8 for 
R|  > 2 ' 10^.  no  separate  trends  could  be  discerned  between  the  two  sets  of  data.  While 
existence  of  differences  may  indicate  some  defects  in  the  analysis  method,  it  is  felt  that  they 
are  small  enough  to  lx-  ignored  when  selecting  the  average  values  of 

Average  values  ot  for  each  model  in  the  absence  of  significant  wavemaking  are  listed 
in  fable  6.  Included  are  values  of  as  defined  previously,  based  on  model  surface  area, 
and  values  of  C^,  based  on  volume  to  the  two-thirds  power,  ^ Ry^,  “ I^r  fjJ  ( Volume 

fach  value  was  based  on  the  faired  curve  - dotted  line  in  figure  13  - obtained  from  all  the 
data  tm  the  appropriate  model  within  the  R|  range  lO’  < Rj  < 2 • 10^.  For  Reynolds 
numbers  less  than  R,  = 10^,  the  experimental  scatter  is  too  large  for  the  data  to  be  useful, 
for  R|  > 2 • lO’,  wavemaking  resistance  becomes  significant. 

TABLE  6 - VALUES  OF  RESIDUAL  RESISTANCE  COEFFICIENTS 


Model 

L,/D 

T 

Cr  • 

10^<Rj^ 

10‘^ 

<2-10' 

S 

Vol 

4620-1 

0.500 

0.667 

0.750 

0,25 

2.1 

4620-2 

1.000 

0.667 

1.500 

0,24 

2.0 

4620  3 

1.820 

0.667 

2.730 

0.20 

1.7 

4620-4 

3.000 

0.667 

4.500 

0.20 

1.7 

4620  5 

1.000 

0.850 

1.176 

0.23 

1.9 

4620  6 

1.820 

0.850 

2.141 

0.20 

1.6 

4620  7 

0.500 

0.850 

0.588 

0.28 

2.3 

4620  8 

1.216 

0.823 

1.478 

0.21 

1.7 

4620  9 

1.000 

0.933 

1.072 

0,26 

2.2 

I able  6 shows  that  the  measured  values  of  Cj^,  based  on  S,  for  all  forehodies  ranged 
between  0.00020  and  0.00028.  Models  3,  4,  and  6 had  the  minimum  ( values  of  0.00020, 
indicating  no  gain  or  loss  in  Cj^  when  increasing  L,  /I)  from  1.82  to  3.00.  For  parent  Model 
3(L,  ,'!)=  1 .82),  increasing  Cp  from  0.6O7  to  0.850  did  not  change  the  value  of  Cj^.  Model 
8 is  a llat-faced  body  with  a smaller  L,  /I)  (1.216)  and  a larger  Cp^,  than  Model  3.  these 
changes  caused  only  a slight  increase  in  the  value  C'j^  from  0.00020  to  0.00021. 
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IX-crciisinp  the  value  of  L|.  D to  1.00  had  the  effect  of  increasing  The  three  fore- 
bodies  tested  witli  this  L,.  /I)  (Models  2,  5.  and  0)  had  Cp  values  of  0.667,0.850.  and  0.933. 
Model  9 with  the  largest  C p value  has  the  largest  value  of  ('  (0.00026).  Model  2 
(C'p^  = 0.667)  had  a C'^  value  of  0.00024  and  Model  4620-5  (C'p^  = 0.850)  had  a value  of 
0.00023.  Decreasing  the  value  of  L,  /D  to  0.50  further  increased  the  values  of  Two 
models,  1 and  7.  had  Lj  D’s  of  0.50  and  C’p^  values  of  0.667  and  0.850,  respectively.  The 
measured  values  of  C'j^  were  0.00025  for  Model  4620-1  and  0.00028  for  Mode*  4620-5. 

The  values  of  ('j^  are  given  in  Figure  1 5 as  a function  of  a forebody  slenderness 
parameter  T = ( Lj  D)/C’p^ . f-or  values  of  T greater  than  2.0  C'|^  assumes  a constant  minimum 
value  of  0.(K)020.  As  T becomes  lower  than  2.0  the  values  of  rise  steadily.  Values  of  T 
for  each  model  are  given  in  Table  6. 

The  values  of  ( based  on  model  volume  to  the  two-»hirds  power  are  also  given  in 
Table  6.  The  volume  C^’s  change  the  relative  merits  of  the  models  somewhat  but  not 
enough  to  be  significant. 

To  further  illustrate  the  significance  of  C'|^  results,  prototype  values  of  = (Dj^ 

+ I)|  (Rj  )) '(  1 2p  S).  and  the  effective  horsepower  ehp  for  nine  unappended  shapes, 
which  are  geometrically  similar  to  the  models,  have  been  calculated, 
where  Dj^  is  residual  drag 

D|  is  friction  drag,  at  Reynold^  number  Rj 
R|  is  Reynolds  number.  l.lv 
i.  is  prototype  length 
U,,  is  prototype  speed 
S is  prototype  surlace  area. 

Figure  l6  gives  the  values  of  ('i^  as  a function  of  Rj . I he  ehp  values  are  given  in  Table  7 
for  six  representative  speeds.  Also  given  in  fable  7 are  ehp  ratios  when  ehp  of  Model  3 has 
bc-en  used  as  the  b;ise.  Strictly  speaking,  this  ratio  is  a function  of  Rj . However,  the  variation 
is  not  discernable  in  the  three  significant  figures  given,  fhese  computations  were  made  using 
the  Schoenherr  friction  line,  and  no  roughness  allowance  was  used,  fhe  usual  assumption 
was  made  of  setting  C’^  equal  f«)r  both  the  model  and  prototype  scales.  A linear  ratio  of 
prototype-to-model  of  15.79  was  used,  resulting  in  a full-scale  length  of  1 10  m for  sha|H‘  3. 
Results  indicate  a maximum  variation  in  ehp  of  4 percent  for  the  nine  shapes  investigated. 
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TABLE  7 - EFFECTIVE  HORSEPOWER  REQUIRED 
FOR  THE  NINE  SHAPES 


m/iec 

Shape 

5.15 

7.72 

10.30 

12.87 

15.44 

18.02 

EHP 

(EHP)3 

EHP  • 10'^ 

EHP  • 10'^ 

EHP  • 10'^ 

EHP  • 10'^ 

EHP  • 10*'* 

EHP  • 10** 

1 

5.44 

1.76 

4.04 

7.70 

1.31 

2.04 

1.02 

2 

5.42 

1.75 

4.02 

7.66 

1.30 

2.03 

1.02 

3 

5.33 

1.72 

3.94 

7.52 

1.27 

2.00 

1.00 

4 

5.37 

1.73 

3.98 

7.58 

1.28 

2.01 

1.01 

5 

5.38 

1.74 

3.99 

7.60 

1.29 

2.01 

1.01 

6 

5.29 

1.70 

3.91 

7.46 

1.26 

1.97 

0.99 

7 

5.53 

1.78 

4.10 

7.83 

1.34 

2.07 

1.04 

8 

5.34 

1.72 

3.96 

7.54 

1.28 

2.00 

1.00 



5.46 

1.76 

4.05 

7.73 

1.31 

2.03 

1.03 

Note:  Linear  ratio  = 15.79. 

CONCLUSIONS 

riie  following  conclusions  have  been  drawn  from  the  results  obtained. 

1.  Foi  the  range  of  conditions  covered  in  these  experiments,  the  occurrence  and 
location  of  laminar  separation  is  accurately  predicted  by  the  Curle-Skan  modified  Thwaites 
criterion.  Separation  insures  the  onset  of  turbulent  How  a short  distance  after  the  separation 
position. 

2.  On  the  four  models  that  did  not  exhibit  laminar  separation,  the  results  did  not  show 
a single  unique  relation  between  the  measured  flow  properties  in  the  transition  regions  and 
the  corresponding  computed  spatial  amplification  ratios  obtained  by  linear  stability  theory. 
On  Model  3.  the  range  of  amplification  factors  that  correlates  with  the  onset  of  transition 
lies  between  e*  and  e’^.  The  corresponding  range  for  Model  4 is  from  e’  to  e*’  and  on 
Model  6 from  e’^  to  e'^.  The  data  taken  for  Model  8 indicate  that  a nearly  constant  value 
of  e’  correlates  well  with  data  at  the  onset  of  transition. 

3.  The  residual  resistance  coefficient  of  the  nine  models  correlates  reasonably  well  with 
the  forel)ody  slenderness  parameter  (T  = (L,. /D)/(('p^  ).  The  results  indicate  that  for  a given 
stern  and  constant  total  volume,  an  increase  in  forebody  fullness  results  ii.  a small  increase 

in  residual  resistance  coeffleient,  once  the  slenderness  parameter  drops  lower  than  about  2.0. 
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The  maximum  spread  in  values  of  effective  horsepower  is  about  4 percent.  Models  3.  4.  6, 
and  8 showed  the  best  resistance  characteristics,  while  Model  7 had  the  most  unlavorable. 
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Figure  I ~ Body  of  Revolution 


Figure  2 — Profiles  of  Forebodies  of  Series  I 
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Figure  3 - Profiles  of  Forebodies  of  Series  2 
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Figure  5 - Computed  Distributions  of  Pressure  Coefficient 
on  AxLsymmetric  Forebodies  of  Series  2 
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Figure  7 - Representative  Hot  Film  Signals 

for  various  Flow  Regimes 
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Figure  K - Idealized  Body  Frictional  Drag  Distribution 
with  a Turbulence  Trip 


Figure  10  - Measured  Flow  Regimes  on  Model  4, 
Compared  with  Amplification  Factor 
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Figure  1 1 — Measured  Flow  Regimes  on  Model  6, 
Compared  with  Amplification  Factor 
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2 NAVSHIPYO  MARE 

1 Library 
1 Code  250 

1 NAVSHIPYO  PUGET/Lib 

1 NAVSHIPYO  PEAR L/Code  202.32 

8 NAVSEC 

1 SEC6034B 
1 SEC  61 10 
1 SEC6114 
1 SEC  6120 
1 SEC  6136 
1 SEC6140B 
1 SEC  6144 
1 SEC  6148 

1 NAVSEC.  NORVA/6660.03  Blount 

12  OOC 

1 AFOSR/NAM 

1 AFFOL/FYS.J.  Olsen 

2 MARAO 

1 Oiv  of  Sliip  R&D 
1 Lib 


1 NAOA  HQS/Lib 
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3 


NBS 
1 Lib 

1 P S.  KlebanoM 
1 G.  Kulin 


fRBBEOXMO  MOI  BUMMK/t  Tl 


Copies 


Copies 


1 NSF/Eng  Lib 

1 LC/Sci  & Tech 

1 DOT/Lib  TAD  491,1 

2 MMA 

1 National  Maritime  Research  Center 
1 Litxary 

1 U.  of  Biiftgeport/E.  Uram 

4 U.  of  Cal/Dept  Naval  Arch,  Berkeley 

1 Library 
1 W.  Webster 
1 J.  Paulling 
1 J.  Wehausen 

2 U.  of  Cal,  San  Diego 

1 A T.  Ellis 
1 Scripps  Inst  Lib 

5 CIT 

1 Aero  Lib 
1 T.V.  Wu 
1 A.J.  Acosta 
1 I Sabersky 
1 D.  Coles 

2 Catholic  U.  of  Amer. /Civil  & Mech  Eng 
1 M.  Casarella 
1 Library 

1 Colorado  State  U./Eng  Res  Cen 

1 U.  of  Connecticut /Scottron 

1 Cornell  U./Shen 

2 Florida  Atlantic  U. 

1 Tech  Lib 
1 S Dunne 

2 Harvard  U. 

I G.  Carrier 
1 Gordon  McKay  Lib 

1 U.  of  Hawaii/Bretschneider 

1 U.  of  Illinois/J.  Robertson 

4 U.  of  Iowa 

1 Library 
1 L.  Landweber 
1 J.  Kennedy 
1 V.C.  Patel 

1 Johns  Hopkins  U. /Phillips 

1 Kansas  State  U. /Nesmith 


1 

1 

5 


4 


4 


2 


2 


4 


1 

2 

5 


1 

1 

2 


U.  of  Kansas/Civil  Eng  Lib 

Lehigh  U. /Fritz  Eng  Lab  Lib 

MIT 

1 Library 
1 P.  Leehey 
1 P.  Mandel 
1 M.  Abkowitz 
1 J.N.  Newman 

U.  of  Minn/St.  Anthony  Falls 
1 Silberman 
1 Lib 
1 Song 
1 R.  Arndt 

U.  of  Mich/NAME 
1 Library 
1 F.  Ogilivie 
1 Hammitt 
1 Cough 

U.  of  Notre  Dame 
1 Eng  Lib 
1 Strandhagen 

New  York  U./Courant  Inst 
1 A.  Peters 
1 J,  Stoker 

Penn  State 

1 B.R.  Parkin 
1 R.E.  Henderson 
1 J L.  Lumley 
1 ARL  Lib 

Princeton  U./Mellor 

U.  of  Rhode  Island 
1 F.M.  White 
1 T.  Kowalski 

SIT 

1 Library 
1 Breslin 
1 Savitsky 
1 P.W.  Brown 
1 Tsakonas 

U.  of  Texas/Arl  Lib 

Utah  State  U./Jeppson 

Southwest  Res  Inst 
1 Applied  Mech  Rev 
1 Abramson 
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Copies 


Copies 


3 Stanford  U 

1  Eng  Lib 

1 R Street,  Dept  Civil  Eng 
1 S.J.  Kline,  Dept  Mech  Eng 

1 Stanford  Res  Inst/Lib 

1 U.  of  Washington/Arl  Tech  Lib 

3 VPI 

1 H.L.  Moses,  Dept  Mech  Eng 
1 O.P.  Telionis,  Dept  Mech  Eng 
1 J.  Schetz,  Dept  Aero  & Ocean  Eng 

3 Webb  Inst 

1 Library 
1 Lewis 
1 Ward 

1 Woods  Hole/Ocean  Eng 

1 Worchester  Pl/Tech  Lib 

1 SNAME/Tech  Lib 

1 Bethlehem  Steel /Sparrows  Point 

1 Bethlehem  Steel/New  York/Lib 

1 Bolt,  Beianek  & Newman/Lib 

1 Dynamic  Technology  Corp 

1 Exxon,  NY/Design  Div,  Tank  Dept 

1 Exxon  Math  & System,  Inc. 

1 Flow  Research  Corp 

1 General  Dynamics,  EB/Boatwright 

1 Gibbs  & Cox/Tech  Info 

4 Hydronautics 

1 Library 
1 E.  Miller 
1 V Johnson 
1 C.C.  Hsu 

1 Lockheed,  Sunnyvale/Waid 

2 McDonnell  Douglas,  Long  Beach 

1 T.  Cebeci 
1 J.  Hess 

1 Newport  News  Shipbuilding/Lib 

1 Nielsen  Eng  & Res 

1 Northrup/Wortman 

1 Oceanics 


3 Rand  Corp 

1 E R Van  Drie't 
1 C Gazley 
1 J Aroesty 

2 Rockwell  International  Autonetics 

1 R Scotti 
1 B.  Carmichael 

1 Sperry  Rand/Tech  Lib 

1 Son  Shipbuilding/Chief  Naval  Arch 

1 Robert  Taggart 

1 Tracor 

2 Westinghouse  Electric  (ANNAP) 

1 M.S.  Macovsky 
1 Gulino 


/ 
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CENTER  DISTRIBUTION 


Copies 


Code 

Copies 

Code 

012 

R.C.  Allen 

1 

1552 

H.T.  Wang 

11 

W.M.  Ellsworth 

1 

1552 

J.  McCarthy 

117 

R.M.  Stevens 

20 

1552 

J.  Power 

1500 

W.E.  Cummins 

1 

1560 

G.  Hagen 

1504 

V.J.  Monacella 

1 

1562 

M.  Martin 

1506 

M.K,  Ochi 

1 

1562 

C M.  Lee 

1507 

D.  CieslowksI 

1 

1564 

J.  Feldman 

1508 

F,  Peterson 

1 

1568 

G.  Cox 

1512 

J.B.  Hadler 

1 

1572 

M.D.  Ochi 

1520 

R.  Wermter 

1 

1572 

E.  Zarnick 

1521 

P.  Pien 

1 

1576 

W.E.  Smith 

1524 

G.  West 

1 

1615 

R.J.  Furey 

1524 

W.C.  Lin 

1 

1802.2 

F.  Frenkiel 

1524 

Day 

1 

183 

E.  Cuthill 

1524 

Scragg 

1 

184 

H.  Lugt 

1532 

G.  Dobay 

1 

1843 

J.  Schot 

1532 

M.  Wilson 

1 

1843 

C.  Dawson 

1540 

W.B.  Morgan 

1 

19 

M.M.  Sevik 

1541 

Granville 

1 

1942 

J.T.  Shen 

1542 

Yim 

1 

1942 

W.R.  Brown 

1544 

Gumming 

1 

1942 

F.C.  DeMetz 

1544 

Boswell 

1 

1942 

B.  Bowers 

1544 

Caster 

1 

1942 

T.C.  Mathews 

1544 

Cox 

1 

1946 

1552 

M.  Chang 

1552 

T.T.  Huang 

30 

5214.1 

Reports  Distribution 

1552 

N.  Santelli 

1 

522.1 

Unclassified  Library 

1552 

C.  von  Kerczek 

1 

522.2 

Unclassified  Library 
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